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The curvature correction term to the surface tension is used as a criterion for the efficiency of 
the Skyrme interaction in describing surface properties. Based on the nuclear equation of state, 
the curvature correction term to the surface tension coefficient is calculated for 97 standard Skyrme 
interaction parameter sets in the vicinity of nuclear saturation density at zero temperature. The 
main idea is to find those parametrizations that give the Tolman’s S correction close to the available 
theoretical predictions from the statistical theory. Only 59 out of 97 models give satisfactory results. 

Comparison of the obtained results with the results of the implementation of different macroscopic 
and microscopic constraints to Skyrme parametrizations available in the literature allows us to select 
4 models that satisfy all the constraints. 

PACS numbers: 21.30.Fe, 21.65.-f, 68.03.Cd 


I. INTRODUCTION. 

The derivation of the equation of state (EOS) of nu¬ 
clear matter is among the most important goals and long¬ 
standing unsolved problems in nuclear physics and as¬ 
trophysics [I|. Various approaches to the description of 
infinite nuclear matter exist. Among them are purely 
microscopic ones based on the realistic description of the 
nucleon-nucleon (NN) interaction in which case the 
result depends not only on the chosen interaction but also 
on the way many-body effects are treated. These may be 
handled either by a direct description of the tree-particle 
interactions or by approaches like that of the Bruckner- 
Hartree-Fock method, the Dirac-Bruckner-Hartree-Fock 
formalism [^, or self-consistent Green’s functions Q , 
etc. 

At the same time, in describing the experimental data 
and in computer simulation, the most widely used models 
are those based on effective density-dependent NN and 
NNN interactions rather than on realistic ones (e.g., the 
models introduced by Skyrme i and Gogny 0). The 
main problem of such approaches is in the infinite number 
of possible sets of model parameters providing satisfac¬ 
tory description of the ground-state properties of stable 
nuclei. 

A large number of various Skyrme-force parametriza¬ 
tions and theoretical models attempting to describe nu¬ 
clear matter and finite nuclei in a wide range of external 
parameters exists. They all were constructed under spe- 


* k.cherevko@univ.kiev.ua 
^ jenk@bitp.kiev.ua 
^ fszhang@bnu.edu.cn 


cific assumptions that reduce their predictive power Q. 
Combined with the indirect mo del-dependent experimen¬ 
tal methods used to evaluate nuclear matter properties, 
it makes the selection of the realistic sets of parameters 
quite a difficult task. 

Presently over 200 sets of Skyrme parameters are 
known from the literature. They result from the anal¬ 
ysis of various observables, leading to different predic¬ 
tions concerning the behavior of nuclear matter away 
from equilibrium. Recently, a number of interesting and 
important papers systematically checking the sets of m- 
rameters for nuclear matter constraints appeared [1, 0 . 
Important work on the construction of new parametriza¬ 
tions with systematic variation of the parameters to im¬ 
prove the precision of the results for some crucial nuclear 
matter properties has been done in Refs. [lol - [T^ . Such 
investigations may result in the improvement of the equa¬ 
tion of state of nuclear matter, applicable in a wide range 
of parameters. Thus, the search for mo del-independent 
constraints connected with the specific properties of the 
nuclear matter is timely and important. 

Such constraints may result from the properties of the 
interfaces. To start with, let us mention that the Droplet 
model of nuclei [l^ plays a special role among macro¬ 
scopic models. It makes possible the description of aver¬ 
age properties of a saturated system, such as a nucleus, 
consisting of two components (neutrons and protons), 
with account for the boundary effects and the presence 
of a diffuse layer. The surface energy and the properties 
of the surfaces in nuclear matter have been studied in 
a number of papers dm. Although the dependence 
of the surface tension (and surface energy) on the sur¬ 
face curvature as well as its impact on different physical 
properties were also studied by several groups of authors 



2 


[illil- Still, for decades it remains one of the most con¬ 
troversial issues in mesoscopic thermodynamics [l^-[ 2 lj . 

In studies of surface properties of nuclei with mass 
number A, the account for the curvature effects is im¬ 
portant. Within the Droplet model this requires the in¬ 
clusion of additional terms proportional to A 3 in any 
expansion concerning the nuclear properties in terms of 
the fundamental dimensionless ratio, given by ^ 3 , 

which is the ratio of the interparticle spacing rg to the 
nuclear radius R [l^. In these studies, effects connected 
with the surface curvature are specified by the curvature 
correction coefficient 03 accompanying terms of the order 
As. In statistical mechanics and in Gibbs-Tolman’s (G- 
B) thermodynamics, of interfaces it corresponds to to the 
Tolman length, called also 6— correction [^. The basic 
parameter d was first introduced by Tolman It is 

equal to the distance between the equimolar surface Rem 
and the surface of tension R at the interphase boundary 

6 = Rem- R- (1) 

According to the G-T theory, the surface tension a of the 
curved interface, in the leading order approximation, can 
be defined as 

a{R) = a^(^l-^ + -->j, ( 2 ) 


from the literature (see Appendix]^. In our analysis we 
included the most popular parametrizations. From large 
families (e.g. BSk,SkSC), we selected several represen¬ 
tative members. The Tolman length was calculated for 
the whole families in which S coincides with or is close 
to that calculated statistically. The results are compared 
with the theoretical predictions. This analysis is inter¬ 
esting as a test of various models regarding their ability 
to describe interphase interfaces. 

Various approaches for the evaluation of the curvature 
correction exist. A method to calculate the Tolman S- 
correction from the EOS of nuclear mater was introduced 
earlier in [ 2 ^. 

In that model, one gets for S from the EOS of sym¬ 
metric nuclear matter with isospin-independent effective 
mass (see Appendix iBl) in the case of T = 0, at normal 
density po'- 


3 po^ 
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where 


W = 
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where R is the droplet radius (equal to the radius of 
the surface of tension [13, HI) and (Too is the surface 
tension of the planar interface. Originally introduced for 
ordinary liquids, it can be defined for any system with 
curved interface of a non-negligible boundary layer 
such as nuclei and nuclear systems with a finite diffuse 
layer [ 2 ^. 

First theoretical estimates of the correction term were 
done by Tolman It appeared to be close to the av¬ 
erage interparticle distance rg. Namely, 5 ~ 0.3 — 0.6rg , 
that, for a nuclear systems, is rg ^ 0.7 fm at normal den¬ 
sity p ^ 0.17 fm“^. Present calculations from statistical 
mechanics yield Tolman length of the order of the inter¬ 
particle distance S tq = 1.14fm [l^. Thus, mathemat¬ 
ically the term ^ in Eq. ([5]) becomes important for the 
systems with i? < 14 fm (^ > 0.1 in m and even more so 
for nuclear systems with R « 0.7A^ = 0.7(277)3 =< 4.6 
fm for heavy nuclei. 

In view of the importance of the curvature correction 
for nuclear systems, checking different sets of standard 
Skyrme parameters as for their ability to reproduce the 
theoretically predicted values for ^-correction becomes an 
important task. Therefore, following the idea of Q in this 
work we present an attempt to use Tolman correction as 
a constraint for different sets of Skyrme parameters. 

II. THEORETICAL MODEL. 

In the present paper the Tolman ^-corrections are cal¬ 
culated for 97 different sets of Skyrme parameters known 


We use the above results to calculate corrections to 97 
sets of Skyrme parameters, using them in testing various 
parametrizations describing surface effects prop erly . Let 
us recall that, following Tolman’s first estimates [^, sub¬ 
sequent calculations based on statistical mechanicsthat 
use minimum information regarding the details of the in¬ 
teraction about the studied liquids but operating with 
dimensionless parameters such as the size of the diffuse 
layer divided by the interparticle distance (e.g. penetra¬ 
ble sphere model in the mean field approximation)yield 
I (5 I in a range from ^rg to |rg (' [^[291 - 1^ and references 
therein). Even though the sign of the curvature correc¬ 
tions reported in the literature varies in different calcula¬ 
tions, the absolute value is approximately the same, vary¬ 
ing within the distance where the density profile faces 
rapid changes. The density profiles in the diffuse layer in 
the dimensionless coordinates are very similar for ordi¬ 
nary liquids and nuclear matter. This allows us to extend 
the results of the ordinary liquids physics to nuclear mat¬ 
ter. Thus, with account for the uncertainty of the theo¬ 
retical values of S, the admissible range for the curvature 
corrections adopted in the work is 

I S 1= (0.3 - 0.6)ro = (0.34 - 0.68)fm. (5) 

In all calculations the surface tension of the semi-infinite 
matter at T = 0 is calculated for the symmetric case with 
Pn = Pp = \p (without Coulomb interaction) within the 
restricted extended Thomas-Fermi (ETF) approach [^ . 
Terms up to fourth order are considered. All the calcula¬ 
tions include effective mass and spin-orbit contribution. 
The function p{z) used to minimize the surface energy is 
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the one parameter Fermi function; 


Piz) 


Poo 

1 + exp{az) 


( 6 ) 


Our calculations yield values for CToo in the range 0.93 
1.21MeV-fm~‘^. The results of our calculation are pre¬ 
sented in Fig. [T] . 

In Table U the parametrization yielding values of 
S within the admissible range m and, at the same 
time, satisfy the constraints of Ref. @ are shown, ap¬ 
pended by the macroscopic properties of the correspond¬ 
ing parametrizations. The label ’’part” indicates that 
the parametrization passed the test of Ref. @ for all but 
one of the applied constraints, and the failure for that 
one constraint is less than 5 per cent. 


TABLE I. Tolman 5-correction and macroscopic properties 
for various Skyrme parametrizations. 5 is in fm; po is fm~®; 
K is in MeV; Co’’ is in MeV fm®. Forces below ’’Extended” 
pass the 5 - correction test with the extended range up to 
interparticle distance 


Skyrme force 

5 

Po 

K Eo Co'’ 

[9] results 

Ska25s20 

-0.74 

0.1746 

210.78 -15.32 -250.23 

-h 

SV-min 

-0.74 

0.1746 

210.78 -15.32 -250.23 

+ 

Ska35s20 

-0.78 

0.16 

230 -16 -272.70 

part 

SkTl 

-0.58 

0.162 

201.95 -15.81 -240.38 

part 

SkT2 

-0.58 

0.162 

201.95 -15.81 -240.38 

part 

SkT3 

-0.58 

0.162 

201.95 -15.81 -240.38 

part 

Skxs20 

Extended 

-0.58 

0.1595 

234 -15.94 -253.67 

part 

LNS 

-0.74 

0.1746 

210.78 -15.32 -250.23 

+ 

SQMC700 

-0.74 

0.1746 

210.78 -15.32 -250.23 

+ 

MSLO 

-0.58 

0.1595 

234 -15.94 -253.67 

part 

SKRA 

-0.58 

0.1595 

234 -15.94 -253.67 

part 

KDEOvl 

-0.58 

0.1595 

234 -15.94 -253.67 

part 


III. RESULTS AND DISCUSSION. 

Only 59 parametrization of those analyzed satisfy the 
constraint on the 5—corrections imposed in Eq. ©■ 
Among them Ska25s20 and SV-minOsatisfy also all the 
criteria of Ref. a, while Ska35s20,SkTl, SkT2,SkT3 
and Skxs20 satisfy all but one constraints of Ref. a^ If 
one increases the admissible range of 5 up to the interpar¬ 
ticle distance, then two more parametrizations satisfying 
all the constraints of Ref. a three more satisfying 
all but one will pass the 5-correction test (see TableU). At 
the same time, as seen in Fig. [l] many parametrizations 
yield values of 5 correction close to the range allowed 
by ®. An interesting but not surprising observation is 
that, while some of the parametrizations of the family 
SkT pass the 5 constraint test, others do not, although 
the parametrizations are based on the same inputs and 
use the same method. This observation may suggest that 
some of the parametrizations pass the test just by chance. 
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FIG. 1. Color online. Tolman’s 5-length. Solid lines show 
the boundaries defined by ®; dashed line correspond to the 
extended range up to interparticle distance 
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The main features of the parametrizations that passed 
the test are shown in (Tab. [TTl Appendix O. As seen 
from the Table, the families with acceptable values of 
the Tolman correction use finite nuclei properties con¬ 
nected with the nuclear surface {e.g. surface properties 
of selected magic and semimagic nuclei, surface thick¬ 
ness, neutron rms radius) as input data. At the same 
time, families that did not pass the test use different in¬ 
put data {e.g. v070 or skz2). It should be also noted 
that parametrizations elaborated for neutron matter {e.g. 
the SLy family) fail to produce an acceptable curvature 
correction to the surface tension in symmetric nuclear 
matter. 

The obtained values appear negative for all chosen 
parametrizations except for SV (see Fig. [!}, which 
means that the surface of tension is located closer to the 
liquid phase with respect to the equimolar surface. 

An important observation is that all parametrizations 
constructed by systematic variations of the parameters 
(e.g. SV and SkS families) do pass the J-correction test. 

Attempting to find an apparent pattern of the forces 
performance in describing properties of the nuclear sur¬ 
faces, we analyzed the dependence of <5 on various force 
parameters and the macroscopic properties characteris¬ 
tics for the tested forces. We also calculated the coupling 
constants found in recent years and aimed to improve the 
relevance of the Skyrme forces to different physical prop¬ 
erties Q and given as linear combinations of individual 
parameters 

= 1^0 + S^3P0“ /yN 

C/ = -jto Q + xo) - (i -h Xs) po“- 

Other combinations of the individual parameters are 
those providing the most compact formulation of the en¬ 
ergy functional and the residual interaction [s^ 


bo = to (1 + \xo) , 

h = (l -h \xo) . 


( 8 ) 


Unfortunately, no particular dependence of the value 
of 5-correction in almost all of the force parameters, 
coupling constants or macroscopic properties was found. 
The only observable correlations are in the slight increase 
of the absolute value of 5 with increasing ti (Fig. [2]), de¬ 
creasing effective mass (Fig. [3|) and increasing absolute 
value of the coefficient (Fig. S]). From standard 
statistical analysis of the forces that pass the curvature 
correction test one may suggest that with probability 0.95 
the most probable values of the above parameters should 
be 


= -245.1 ±3.6 

= 0.96 ± 0.02 (9) 

ti = 304.1 ±10.2 

The standard deviations are SDcgP = 13.8, SD^* = 0.07 
and SDt^ = 36.7. One can see from the pictures ©-(HI) 
that there exist a dependence of 6 on ti, m* and Cq^. 
Even though this may be the necessary condition for the 
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FIG. 2. Dependence of absolute value of 5 on ti coefficient of 
the Skyrme force 
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FIG. 3. Dependence of absolute value of S on effective mass 
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FIG. 4. Dependence of absolute value of 5 on Co^ 
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force ability to describe nuclear surfaces, certainly it is 
not sufficient, since some of the parametrizations, includ¬ 
ing those with S outside our range ([5]), produce ti, m* 
and Cq^ within the ranges imposed by Eq. This may 
be an indication of some basic problem in the Skyrme- 
type parametrizations, maybe connected with the large 
freedom in choosing the force parameters. 


IV. CONCLUSIONS. 


In this paper we calculated the curvature correction 
term of the surface tension from the nuclear equation of 
state for 97 different standard Skyrme parametrizations 
available in the literature. The obtained results show 
strong dependence of the curvature correction term on 
the EOS. 

To summarize, our study shows that not all the exist¬ 
ing parametrizations are capable of describing adequately 
the interphase interfaces in nuclear matter. In spite of the 
considerable uncertainty regarding the absolute value of 
S, the suggested constraint, even with a wide admissible 
range for 6, allows one to test various Skyrme forces re¬ 
garding their capacity to describe curvature effects. It 
should be mentioned that better agreement is observed 
with the use of EOS that account for the surface effects 
with respect to those that do not (or were designed for 
neutron matter). 

Comparison of the obtained results with the avail¬ 
able data on different Skyrme forces with nuclear mat¬ 
ter constraints suggests two parametrization, namely 
Ska25s20 and SV-min, that satisfy all constraints. For 
an extended admissible range for 6 that number in¬ 
creases and two more parametrizations, namely LNS 
and SCMC700, come into play. It can be seen from our 
study that systematic variation of the parameters is quite 
efficient, giving promising results in studies of the basic 
properties of nuclear matter. At the same time, further 
progress with Skyrme type forces requires a better under¬ 
standing of the physical meaning of different parameters 
rather than the introduction of new parameters. 
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Appendix A: Parametrizations used 


The Skyrme parametrizations analyzed in the work are 
BSkl [si, BSklO [13, E ra, Ecr [T3 i, F^on [sl, Gcr 
KDEOv [H^DEOvl 
MSkS [li^Skd 
MSkASi^SLO 14 SIRAPR 



^1, PRC45 



_J, RATP [il^cr [ii, SGII 
il, SII Hi, SIII 1^, SIV0, SK255 SK272 M, 

Ska [4i^ka25s20 M, Ska35s20 1], SkI4^|, SkM* [5g, 
SkMP li, SkO Isa, SkP Jp, SKRAjp, SkSl 


SG02 


HI 


SkS2 [lJ7skS3 [Tjr SkS4 nr SkSCl HfTskSClO [HI, 
SkSGrH, SkSG3 [53^kSC4 [ll^kSCS [pL SkSC6 
p, SkSPl & SkTl 0, SkTl* tsi, SkT2 |5J, SkT3 
Ml, SkT3* 13, SkT4_|53, SkTSjH, SkTbJsi, SkT7 



Skxs20 
SLv6 


M|, SkT8 
M|, Skxsl5 [& 

SLy4 [ 53 , SLy5 
0 , S(^C700 0, SV^ 
K226 111, SV-kap02 fTT" 
SV-sym32 [p SV-tls [R 

a , ZRlb^, ZRlc^ 

, ZR3a gil, ZR3b gF 


SKX [55|, SKXce 


Skxs25 

, SV-bas E 
SV-mas08 |T 
T [13, v070 [5 


155|, SKXm 
Skz2 [ 53 , 
SQMC650 


ZR2a 

ZR3c 



ZR2b 14 

r Zcr [13 


ZRla 

ZR2c 


Appendix B: Equation of state of nuclear matter 


To calculate the ^correction, an EOS of nuclear mat¬ 
ter at low-temperatures and in the high-densities limit 
was used, where A^p » 1 (i.e. when the average de 
Broglie thermal wavelength A is larger than the average 
interparticle separation p~'^). In this case, the EoS takes 
the form [^ : 


Pip,, T) = i: [le%g{p„ T) - SkM, ^)] 

+ 1 (1 + f+ || (1 + a) (a + l)p«+2 

(xo + 5) - l| (5+2:3) {a + l)p“X;Pg^ 


with 


^kq — 


1 23 



(Bl) 

‘9 P^iVq), 
^P^iVq), 

(B2) 


where m and m* are the mass and effective mass respec¬ 
tively, T and p are temperature and density, q is the par¬ 
ticle type (g^proton, neutron), F is the Fermi integral, 

A = is the average de Broglie thermal wavelength, 

g = 2 is the spin degeneracy factor, to, to, xo, xo and a 
are the Skyrme force parameters and 13 = ^. 


Appendix C: Methods and input data used in 
construction of different Skyrme parametrizations 
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TABLE 11. Comparison of the methods and input data used in construction of different Skyrme parametrizations 


Skyrme forces 

Method 

Input data 

KDEOvl 

Simulated annealing method 

Ground state properties of normal and exotic nuclei: 

binding energy, charge radii and spin-orbit splitting, 

radii for ld 5/2 and I/ 7/2 neutron orbits in *^^0 and "^^Ga 

breathing mode energies, critical density pcr, 

positive slope of the symmetry energy up to 3po, 

enhancement factor associated with GDR, Landau parameter G'q 

LNS 

Brueckner-Hartree-Fock 
with 2- and 3-body forces 
(homogeneous matter) 
Hartree-Fock 
(finite matter) 

Nucleon effective mass in symmetric nuclear matter (SNM) 

and asymmetric nuclear matter (ANM), energy per 

particle in SNM and ANM as function of density and 

proton neutron asymmetry; 

constraint on the Landau parameter Go, 

surface properties of selected magic and semimagic nuclei, 

spin-orbit splitting pi /2 — 4 * 3/2 in 

SV family 

Hartree-Fock 

Properties of finite nuclei (energies, radii and surface thickness), 
energies of giant resonances, 

systematic variations of selected nuclear matter properties 

MSkA 

Density dependent Hartree-Fock 

Empirical binding energies and charge radii of the: 
closed shell nuclei ^®0, "^^Ga, ®°Zr, and ^°®Pb; 
isotopes and ^^'*Pb; doubly closed ^®^Sn; 

SK255, SK272 

Hartree-Fock based random phase Nuclear binding energies, charge radii, and neutron radii; 
approximation (RPA) approach rms charge radius of ^°®Pb; incompressibility coefficient from 

relativistic and nonrelativistic mean field based RPA calculations 

SkS family 

Hartree-Fock seniority 

Symmetric nuclear matter ground state properties; 
binding energies of ^®0, '‘°’‘^®Ga, ®°Zr, and ^°®Pb; 
spin-orbit splitting in ^®0; surface parameter and 
symmetry properties; fission barrier of ^"^^Pu; 
restrictions on Landau parameter for SkS4 

SkT family 

Hartree-Fock; 

Extended-Thomas-Fermi 

Nuclear radii and binding energies; droplet model parameters; 
experimental masses, charge radii, charge distribution, neutron 
skin thickness, semiclassical fission barriers and 

Landau parameters 

Skxs family 

Skyrme Hartree-Fock; 
relativistic mean-field models 

Binding energies, rms charge radii, and single-particle energies; 
binding energy difference ‘^®Ni-‘^®Ga; charge density of ^°®Pb; 
constraint of Q!=1/6 for the density dependent potential 

E, Z, T, Ect, Go-, Ro-, Zct 

Hartree-Fock-Bogoliubov 

Binding energy, diffraction radius, 

surface width of '‘°-‘‘®Ga, ®®Ni, ®°Zr, 208 p^. 

1 — s splitting of the Ip level in 
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